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Abstract 

It has been suggested that the suppression of Jjtp production in heavy nuclei 
is a signature of the formation of quark-gluon plasma. We here show that this 
phenomenon can be understood in terms of conventional physics, i.e. i) perturba- 
tive QCD, a) the parton recombination implementation of shadowing in the initial 
state, and in) final state interactions with the hadronic debris of the nuclear tar- 
get. Unlike previous calculations we include both the direct J/ip production and 
its production via radiative x decays {xj ^ J/i^ + l)- We are able to reproduce 
the experimental data including their small-x behavior. We emphasize the impor- 
tance of studying the X2-dependence of the ratio (T{bA)/a{bN), where b designates 
the beam and X2 is the momentum fraction of the parton from the nuclear target. 



1 Introduction 



It has been suggested that the suppression of J /if) production on heavy nuclei is 
a signal for the formation of the quark-gluon plasma; however we would like to 
consider a conventional physics solution to the problem. J/ip production will be 
described by conventional perturbative QCD after including the modification of 
the small-x behavior of the distribution functions for gluons and sea quarks in 
the nuclear medium [|l|. This allows for the continued use of the hard scattering 
amplitudes and the factorization theorem. This modification, plus a correction for 
the 'classical EMC effect at intermediate x, provides a multiplicative modification 
of the parton distribution functions, for partons in a nuclear target. Furthermore, 
because of the hadronic nature of the final state, effects of the interaction of the 
final state with the hadronic debris of the nuclear target must be included. We 
will show how such a picture can reproduce in detail a wide range of data on lepto- 
and hadroproduction of J/ip. 

J/V', ^' and T suppression on nuclear targets H provides valuable data for 
the study of the A-dependence at small-x region and/or large xp- Here x is the 
Bjorken-x variable, the momentum fraction of a parton within a hadron, and Xp 
is the Feynman-x variable, xp = Xi — X2 where Xi and X2 are Bjorken-x variables 
for the partons from the beam and target, respectively. The experimental energy 
ranges from 40 to 800 GeV, which spans part of the intermediate x range through 
the small-x region for J/ ip production. 

J/ip production in Deep Inelastic Scattering (DIS) has been measured by the 
NMC collaboration for incident muon energies of 200 and 280 GeV. The results 
are presented as a ratio of cross sections for two different targets: Sn and C. 
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These data are relevant to the study of the ^-dependence of the cross-section at 
small X. In fact, A-dependence in J/ip production is already known for hadronic 
processes 0] and the cross section ratio 

has been studied considering initial and final state nuclear effects ^ . This ratio 
is for J/ ip production less than unity in the small x region. This result cannot be 
explained by QCD parton model, which gives the ratio equal to 1. 

Nuclear effects have been observed in high energy processes at different momen- 
tum transfer, Q^. Those effects can be traced in DIS with neutrinos or charged 
lepton beams p, ||, ^ |T0[, in Drell-Yan processes |]1T|, and in hadroproduction of 



heavy quarks |^, [T^, [T^. Only the intensity of the suppression differs, and this 
is a key aspect. The general behavior is Agff/A < 1, where Agff is defined as 
a{bA)/a{bN), with b designating the beam. Common features include a more pro- 
nounced effect for heavier target nuclei, a rapidly diminishing effect with increasing 
X and very little dependence on Q^. The first result on nuclear dependence were 
obtained by EMC at intermediate x |T^. Since then the experimental results have 



been extended to smaller x values, and they exhibit the shadowing phenomena. 
The available data for J/tp production are presented in Table 1. 

2 The model 



In DIS as well as in hadroproduction, the gluon fusion process is dominant [|T5|, |T6 
and the small-x behavior of the target gluon is crucial. In order to understand the 
heavy meson suppression with nuclear targets, it is critical to understand small- 
X behavior of the gluon structure function. We here consider a recombination 



model to take into account shadowing at the parton level in the initial state of 
the process [|I|. This approach introduces a modification of the parton evolution 
equations in order to take into account the superposition probability when the 
partons have a large longitudinal size (or large This model incorporates the 

recombination through ladder diagrams as a perturbative mechanism enabling a 
factorized calculation for the cross section ratios M. This approach successfully 



explains both EMC and Drell-Yan small-x data |T^ . 

The recombination effect is enhanced in the nuclear medium, where the longi- 
tudinal size of the parton, Az, can exceed the size of the nucleon at small- a; region. 
The quantity measured experimentally is a ratio of structure functions, 

and this quantity is found to be approximately independent. The x- and A- 



dependence of this ratio has been parameterized by Berger and Qiu 1T8[ , and these 
authors find that it factors in the DIS case: 

Remc{x, A) ^ Rg{x, A)Ra{x, A). (3) 

The Rg{x,A) factor is associated with the partonic shadowing, and it has the 
functional form: 



1 Xc < X < 1 

Az - Az, " 
AzA - Az, 



Rg{x,A) = \ l-KgiA^^-l 



Xa < X < Xc (4) 



1- Kg{A^/^ -1) 0<x<xa 
where Kg parameterizes the amount of gluon shadowing, Az = l/{xp) is the 
wavelength of the gluon, Azc = 1/ {xcp) is the longitudinal distance at which 
neighboring nucleons begin to interact and Aza = ^/{xap) is the longitudinal 
size of the nucleus. Thus the variables xa and Xc are related to the Bjorken-x 



corresponding to a probe of the nucleus and nucleon, respectively. We assume 
that xa = Xc /A^/^, following eq.(9) of Ref. ||T8[. The remaining factor, Ra{x,A), 



parameterizes the classical EMC effect, and it has the approximate form: 

i?,(x,A) = - + ir,(l--), (5) 

Xi X\ 

with this parameterization valid for < a; < 0.6. 

The above factors incorporate the 'classical EMC effect and initial state effects. 
Final state effects are incorporated via a factor, Rss{xf,A) = A^^^^^)"^), where 
a^xp) = 0.97 — 0.27x|^. This expression is suggested by data and also agrees with 
open charm results. This should be adequate. It is premature to deal with specific 
final state effects without having a better understanding of the gluon distribution 
function behavior at small-x. Also, our parameterization of final state effects 
depends only on a final state variable (xp)- 

Other descriptions of final state interactions have been suggested. The nuclear 
approach a la Glauber includes rescatterings of the QQ or heavy meson in the 



nuclear medium |T9[. Other authors also consider specific final state effects (and 
attempts to make quantitative predictions of them) |Q or Pomeron exchange 
models PTl, or raise the fundamental question of the validity of factorization in 



this kinematical region ||2^, |23|. 

3 J/i/j in DIS 

Deep inelastic photoproduction of J/ is understood to take place by the photon- 

gluon fusion mechanism []T6[ 7 + (71 J /ip + g2- The extra gluon in the final state 



is required for the cc to be produced as a color singlet with the correct quantum 
numbers [J^ = 1^) of the J/ip; color and spin projection techniques were used to 



extract the relevant part of the 7 + (7 — cc + amphtude. It has been shown that 
in the inelastic region {z = Ej/^/E^ < 0.9) both gluons are hard and it is not very 
relevant to take care of higher order multiple gluon diagrams. 

We combine the color singlet model for leptoproduction (photoproduction) of 
the J/ ip with the parton recombination model to account for the A-dependence of 
the gluon function. We use the Weizsacker- Williams approximation for the photon 



and the Morfin-Tung leading order set of parton distributions |^J] for the gluon. An 
advantage of DIS compared to hadroproduction in that there are no interactions 
with hadronic components of the beam . 



The color singlet model for J/ip production reproduces the rapidity distribu- 
tion The color singlet model combined with the Weizsacker- Williams ap- 
proximation has been shown to agree with the electroproduction of J/tp with a 



15 GeV electron beam at SLAG |26| and muoproduction with a 280 GeV muon 



beam at GERN |2^, as long as experimental cuts ensure that the J/ip produc- 



tion is inelastic. We therefore feel confident in applying the method to our study of 
shadowing in J/tp muoproduction on nuclear targets. Also, the model has the in- 
teresting features of providing a direct measure of the gluon distribution function, 
in the case — > J/ip + X, or the direct determination of the gluon A-dependence 
by means of Eq. (1). 

We assume that there is no EMG effect in the carbon target, since it is light and 
for tin, the parameters Ka and xi are 1.20 and 0.25 respectively. We also choose a 
fixed value of Xc for carbon, given by Eq. (9) in Ref. |]18|, which for A = 12 gives 



Xc = 0.10. For tin, we allow Xc and Kg to be free parameters, and try to fit the 
existing data. We choose the ranges < Xc < 0.25 and < Kg < 0.50. We then 
calculate the ratio of the cross sections at the x values of the NMG data points, 
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and perform a analysis. We now describe our results for A-dependence on DIS 
and hadropro duct ion of J/ip 

The DIS data is still not abundant but it is extremely useful for comparison 
with the hadronic case. Our results are presented in Fig. 1 and it is clear that 
more statistics are needed to enable a more critical analysis of this model. Due to 
the large error bars it is difficult to constrain the parameters. Also an extension to 
lower X is important to provide a better definition of the x behavior of the ratio, 
Eqn. (2). 

4 J/i/j in hp{A) 

The hadroproduction of J/ip can proceed through a number of parton level pro- 
cesses. The leading order source of J/ip in hadron-hadron collisions is due to 



gluon-gluon fusion {g + g J/^p + g) p8|, although it is not the most im- 
protant contribution. The combination of a lower order (in a^) and large branch- 
ing fraction makes xj production, followed by radiative decay (xj J/i^ + l): 
the dominant source of J/ip (hereafter referred to as radiative production). The 
leading contribution to xj production is the low pt process g + g ^ Xo,2- Ad- 
ditional contributions of the same order in as as direct production come from 
g + a ^ Xo,i,2 + g, q + g ^ Xo,2 + g and qq Xo,2 + g- All the required parton 



level cross sections can be found, e.g., in Refs. [ER |2S|. Furthermore, we use the 



Morfin-Tung leading order parton distribution functions p4[ for partons from the 



target or from the proton beam, and Owens pion set 1 for partons from the 
pion beam. We are attempting to perform a more careful calculation of J/ip pro- 
duction on heavy nuclear targets than previous analyses, and so the cross sections 
we use are those in which the correct color singlet structure, the correct angular 



momentum quantum numbers and small relative momenta are projected out.Q We 



then include the factors Remc and Rss in the calculation of the various cross sec- 
tions for production on the heavy nuclear target. The J/ip cross section and X2 
distribution can be constructed from the direct and radiative xj components, with 
the inclusion of branching ratios where appropriate. 

As in the DIS case, we allow and Kg to be free parameters, and try to fit 
the existing data. We choose the ranges < Xc < 0.25 and < i^^ < 0.50. 
We then compare the experimental data points on R{x2) (where X2 is the parton 
momentum fraction of the gluon from the target nucleus) to our calculations, and 
perform a analysis. 

We consider data from p and vr beams with energies from 200 to 800 GeV, 
focusing on higher mass nuclear targets where the effect of a heavy nucleus is more 
pronounced. In Fig. 2 the results are presented for J/ip (Fig. 2a) and ip' (Fig. 2b). 
The best agreement requires different parameters for each case. This should not 
be surprising since the i})' is a spatially bigger resonance than J A more refined 

version of this calculation should include this fact in the final state effects. We 

^Because of the existence of g + 5 ^ Xo,2, the higher order (non-zero pt) processes involving 
Xo,2 diverge at low pT- This is merely an artifact of an incomplete calculation. If one performs 
a full calculation of xo,2 (5 + 3 ^ Xo,2 including all 1 loop graphs) production, and cancels 
the low Pt divergences against virtual infrared divergences and then absorbs the remaining 
coUinear divergences into the running of the parton distribution functions, the cross section is 
indeed finite. We adopt a less rigorous approach to this problem. It is known that at low pT, 
da /dp\ oc e~^*^^ where Mt is the transverse mass of the charmonium state. Also, the divergences 
in \A(gg Xo,25)P and \A{qg Xo,2'?)P are 1/t, which gives l/p\ at low pr- Therefore we 
regularize these squared amplitudes with a factor 



where ptq is a free parameter. At low pT, \A\'^ should be (a slowly varying function of pt)/Ptj 
and so our regularization should reproduce the observed p'^ distribution. To find a value for pro, 
we first calculate the total xo,2 cross section from the g + g ^ Xo,2 subprocesses and vary pro 
until the integration of the regularized \A\'^ yields the correct value. 




(6) 
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obtain very good agreement for 800 GeV for the J/ case. It is interesting to note 
that for the ip' case, the value of Xc required for the best fit is rather large (the 
largest we allow in our analysis). A reasonable fit to the data can be obtained, 
in this case, in a range of parameter space as demonstrated by the dashed line in 
Fig. 2b. 

Also for 200 GeV and p beam the agreement is good and the results are shown 
in Fig. 3. However for tt we reproduce quite well the general behavior at 280 and 
200 GeV, as shown in Figs. 4 and 5, respectively, but for smaller x our result is 
below the data. It is likely that the pion distribution function used is not suitable 
for this kinematical region. 

As an overall result, considering that higher statistics are still needed to clarify 
this complicated problem, we believe that a conventional model as the one pre- 
sented here is able to accommodate the data. However, it should be emphasized 
that a better understanding of the small-x behavior of the gluon distribution func- 
tion is needed, and in this HERA and photoproduction experiments can play an 
important role. 

5 Discussion 

The inclusion of the finite pt contributions modifies the kinematics of the problem 
somewhat. Previous analyses have assumed that a 2 ^ 1 subprocess dominates 
the J/ip production. If this is the case, it is apparent that a fixed cm energy {^/s) 
and a fixed invariant mass (M^), implies a fixed r = M^/s. But since r = ^1X2, the 
relations between the kinematical variables is given by Xi = {\J^t + + xp)/2^ 
X2 — {\J^T -\- Xp — xf)/2. Simply putting r = Mj^^/s and measuring xp, one 
can extract the parton momentum fractions. Now, however, with the possible 
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addition of more final state partons, the invariant mass of the produced state is no 
longer the J /ip mass, and so the kinematical relations must be modified somewhat. 
In Fig. 6, we show the invariant mass distribution for J/ip and T. The average 
invariant mass in J /if) production is over 1 GeV above the J/ip mass (when we 
include only the g + g J/ip + g subprocess), but the difference is relatively 
smaller in the T case, so use of the correct kinematical expressions for xi and X2 
is in order. The inclusion of radiative xj decays will not significantly alter the 
preceding argument. The correct expressions for xi and X2 depend on the mass of 
the produced charmonium state, it's pt and energy measured in the lab frame, and 



xp- The expression can be simply derived starting from equation (3.2) in Ref. |jT6 

M2 



+ 



z(l-z) 



(7) 



where z = Eonia/Eg^ with the energies measured in the lab frame and gi is the 
gluon from the beam. Replace z with Zobs/xi {zobs = Eonia/Ebeam) and s with 
X1X2S and solve for xi and X2 in terms of experimentally measurable quantities. 



Xl,2 



{xf + z, 



obs ) 



± 



\ 



(xf + ZobsY + 



^^obs 



2[xFM^-Zobs{pl-M^)] 



SZ(^bs 



where Mt = yPr + is the transverse mass of the charmonium state produced. 
It is important to note that the values for xi and X2 assuming fixed r are not even 
correct at pt = 0, since it is possible for the final state parton to be collinear with 
the charmonia state without being soft. 

The small-x region represents one of the last frontiers of perturbative QCD p2| . 
For this reason as well as for the very important study of gluon shadowing more 
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data, both at lower x and with higher statistics, are needed. 
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Exp crime lit 


Beam 
energy 
(GeV) 


Beam 
type 


Target 


2^2 

range 


Xp 

range 


Reference 


AT A O 

NA3 


43 


TT" 


Be, Cu, W 


0.1 - 0.3 


- 0.95 


rFTi" 1 

m 


AT A O 

NA3 


39.5 




H, W 


0.11 - 0.34 


- 0.85 


[ |r> rJl 


AT A O 

NA3 


150 
200 
280 


TT", p 


H, Pt 

TT T~» 1 

H, Pt 
H, Pt 


0.03 - 0.17 
0.0225 - 0.15 
0.016 - 0.127 


0.01 - 0.9 
0.01 - 0.8 
0.01 - 0.9 


fin 

[13 


E537 


125 




TUT T ~) 

W, Cu, Be 


0.04 - 0.20 


0.02 - 0.75 


i 


E672 


530 


TT^ 


C, Al, 
Cu, Pb 


0.013 - 0.124 


0.1 - 0.8 


[0 


E772 


800 


P 


D, C, Ca, 
re, W 


0.01 - 0.4 


0.15 - 0.65 


i 


17770 


snn 

oUU 


P 


u, oa, 
Fe, W 


U. i - u.oo 


n 7 
u - U. ( 


iuM 
(T) 


F7nfi 

Hi i \J0 


OUU 


P, P) 




n ni n 1 99 

u.uiu - \J.l^^ 


u - u.^o 




NMC 


280 




H, D 


0.02 - 0.3 


NA 


m 


NMC 
NA37 


200 
280 




C, Sn 


0.02 - 0.2 


NA 


i 



Table 1: Summary of experimental results on J/ip production on nuclear targets. 
The values of X2 were derived using r = Mj^^/s. 
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Figure Captions 



Figure 1 - Production of J /ip for a /i beam of energy 280 GeV on tin {Sn) and 
carbon (C) targets (experimental data from Ref. 0). The curves correspond to 
shadowing {xc = 0.03 and Kg = 0.50), EMC effect and strong screening (dotted), 
shadowing {xc = 0.10 and Kg = 0.20 for comparison), EMC effect and strong 
screening (dashed), shadowing (xc = 0.03 and Kg = 0.50) and strong screening 
(dotdashed) and strong screening alone (sohd). 

Figure 2 - Production of J/ip (a) and ip' (b) for a proton beam of energy 
800 GeV on tungsten {W) and hydrogen (-^2) targets (experimental data from 
Ref. 0). In Figure 2a, the best fit requires the shadowing parameters Xc = 0.185 
and Kg = 0.04, while in Figure 2b, the best fit requires Xc = 0.25 and Kg = 0.10 
(solid hne). Also shown Figure 2b is the fit for Xc = 0.20 and Kg = 0.135 (dashed 
line) . 

Figure 3 - Production of J/ ip for a proton beam of energy 200 GeV on platinum 
{Pt) and if 2 targets (experimental data from Ref. [|13]). The curve shown includes 
only final state effects and the 'classical EMC effect. 

Figure 4 - Production of J /ip for a vr" beam of energy 280 GeV on Pt and H2 



targets (experimental data from Ref. [jl3[). The best fit is achieved by including 



only final state effects and the 'classical EMC effect. 

Figure 5 - Production of J /ip for a vr" beam of energy 200 GeV on Pt and H2 



targets (experimental data from Ref. [13|). The best fit is achieved by including 



only final state effects and the 'classical EMC effect. 

Figure 6 - The cross section vs. invariant mass (V^) for the subprocess g+g 
J/ip + g, assuming a proton beam of energy 800 GeV on a H2 target. 
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